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ABSTRACT 

We have modelled Atacama Large Millimeter/sub-millimeter Array (ALMA) long 
baseline imaging of the strong gravitational lens system H-ATLAS J090311.6+003906 
(SDP.81). We have reconstructed the distribution of band 6 and 7 continuum emission 
in the z = 3.042 source and we have determined its kinematic properties by recon¬ 
structing CO(5-4) and CO(8-7) line emission in bands 4 and 6. The continuum imaging 
reveals a highly non-uniform distribution of dust with clumps on scales of ^ 200 pc. 
In contrast, the CO line emission shows a relatively smooth, disk-like velocity field 
which is well fit by a rotating disk model with an inclination angle of (40 ± 5)° and 
an asymptotic rotation velocity of 320 kms -1 . The inferred dynamical mass within 
1.5 kpc is (3.5 ± 0.5) x 10 10 M© which is comparable to the total molecular gas masses 
of (2.7 ± 0.5) x 10 10 M© and (3.5 ± 0.6) x 10 10 M© from the dust continuum emission 
and CO emission respectively. Our new reconstruction of the lensed HST near-infrared 
emission shows two objects which appear to be interacting, with the rotating disk of 
gas and dust revealed by ALMA distinctly offset from the near-infrared emission. The 
clumpy nature of the dust and a low value of the Toomre parameter of Q ~ 0.3 sug¬ 
gest that the disk is in a state of collapse. We estimate a star formation rate in the 
disk of 470 ± 80M©/yr with an efficiency ~ 65 times greater than typical low-redshift 
galaxies. Our findings add to the growing body of evidence that the most infra-red 
luminous, dust obscured galaxies in the high redshift Universe represent a population 
of merger induced starbursts. 
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1 INTRODUCTION 


Our understanding of high redshift sub-millimetre (submm) 
bright galaxies (SMGs) has grown im mensely sinc e their 
* E-mail: simon.dye@nottingham.ac.uk discovery nearly two decades ago (|Smail et al. I Il997l : 
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iHughes et al. |[l998l : iBarger et al. 1119981 ). The fact that ap¬ 
proximately half of the total energy output from stars within 
the observable history of the Universe has bee n absorbed by 
dust and re-emitted at submm wavelengths (|Puget et al. I 
ll996l : lFixsen et al. Ill998h and that SMGs represent the most 
active sites of dusty star formation at high redshifts, indi¬ 
cates that their role in early galaxy formation is an impor¬ 
tant one. 

Morphological and kinematical measurements of SMGs 
have led many studies to conclude that they are a more en- 


ergetic version of more local ultra-luminous infrared galax- 

ies (ULIRGs; e.g. Engel et al. 

20101: Swinbank et al. 2010: 

Alaghband-Zadeh et al. 2012l: 

Rowlands et al. 201415. How- 


ever, observations have been limited by the low imaging reso¬ 
lution typically offered by submm facilities, forcing detailed 
investigations to turn to other wavelengths which permit 
higher resolution. With strong attenuation over ultravio¬ 
let to near-infrared wavelengths, where imaging technology 
has benefited from a longer period of development, this has 
proven challenging. Hence, a reliance has traditionally been 
made on correlations with other wavelengths which often 
only result in indirect diagnostics of the internal energetics 
of the physical processes at work in these galaxies. 

A powerful diagnostic which gives unique insight 
into the star formation processes in galaxies in gen- 
e ral is measurem ent of molecular gas (see for example 
IPapadopoulos et al. I l2014i . and references therein). In par¬ 
ticular, in the more extreme environments of ULIRG and 
SMG interiors where strong feedback from star-format ion 
and shock-heating of molecular gas by supernovae are dom¬ 
inant processes, star formation mo dels can be p ut th rough 
the most rigorous of tests (e.g. IPapadopoulos et al. I l201lh . 
In this way, examination of the quantity and kinematical 
properties of molecular gas provides a direct probe of the 
mode of star formation. Following this approach, some stud¬ 
ies have concluded that star formation mechanisms in early 
systems are distinctly differ ent fro m tho se see n in mo re lo¬ 
cal sy stems fe.g. lBournaud &; Elmegreen ll2009l : [lo nes et al. I 


A detailed study to measure the properties of star 
formation in distant SMGs requires two key ingredients. 
Firstly, observations must be carried out at submm wave¬ 
lengths, where most of the bolometric luminosity is emit¬ 
ted, to allow emission from the dust enshrouded molecular 
gas to be detected. Secondly, the observations must be of 
sufficient angular resolution to isolate the dynamics of the 
~ 200 pc gravitationally unstable regions u sually found in 
their star-forming disks ([Downes &; Salomonv Ill998h . 

To provide a sample of suitable SMGs for such inves¬ 
tigation, the large area surveys carried out recently by the 
Herschel Space Observatory, such as the Hersch el As trophys- 
ical T erahertz Large Area Survey (H-ATLAS; lEales et al. I 
l201Ch and the Herschel Extragalac tic Multi-tiered Extra- 
galactic Survey ([Oliver et al. I l2012h along with the survey 
at mill imetre wavelengths carr i ed out at the South Pole Tele¬ 
scope ([Carlstrom et al. Il201ll : IVieira et al. Il20l3 i now pro¬ 
vide a bountiful supply of high redshift dusty star bursts for 
detailed study. Obtaining the required high resolution imag¬ 
ing in the submm is now made possible using the Atacama 
Large Millimetre/sub-millimetre Array (ALMA). 

A particularly compelling use of ALMA for these pur¬ 
poses is to target strongly lensed SMGs. The submm 


has long been suspected to harbour a rich seam of 
strongly lensed galaxies due to a high magnif ication 
bias r esulting from t heir steep number counts ain I 
Il996l : INegrello et al. I l2007h . Thanks to the aforemen¬ 
tioned mm/su bmm surveys , such suspicions hav e now 
b een ver ifie d (| Vieira e t al. 1201(3: INegrello et al. I l20ld : 
iHezaveh fc Holder 11201 ll : I Ward low el a.l.lr2013ll . The intrin- 
sic flux and spatial magnifications by factors of ~ 10 — 30 
inherent in strongly lensed systems therefore combine with 
the use of ALMA to provide the highest possible resolution 
and signal-to-noise imaging of SMGs currently achievable by 
a considerable margin. 

In this paper, we report our analysis of the recently re¬ 
leased ALMA science verification observations of the strong 
lens system H-ATLAS J090311.6+003906 (SDP.81), one of 
the first five strongly lensed su bmm sources detected in the 
H-ATLAS data ([Negrello et al. Il2010h . The system was sub¬ 
sequently followed up in th e near-infrared u sing the Hubble 
Space Telescope (HST; see iNegrello et al.ll20l4 for details 
of these o bservations; N14 hereafter). Lens modelling of the 
system bv lDve et al. I (12014 D14 hereafter) showed that the 
observed Einstein ring can be explained by a single compo¬ 
nent of emission in the source plane. The purpose of this 
paper is to exploit the high resolution ALMA imaging of 
the highly magnified lensed source to determine its physical 
properties. One of our key questions is how the rest-frame 
optical source emission reconstructed by D14 relates to the 
reconstructed submm emission detected by ALMA. 

The layout of this paper is as follows: Section [2] out¬ 
lines the data. In Section [3] we describe the modelling pro¬ 
cedure used to obtain the results which are given in Sec¬ 
tion 01 We discuss our findings in Section [5] and summarise 
the major results of this work in Section [6] Throughout 
this paper, we assume the following cosmological param¬ 
eters; Hp = 67 kms~ 1 Mpc " 1 , — 0.32, = 0.68 

([Planck Collaboration ll2014l h 

2 DATA 

2.1 ALMA data 

ALMA Science Verification data on SDP.81 were taken from 
the ALMA Science PortaJ]] (ASP). We give an overview 
of those data here, although more d etails can be found in 
I ALMA Partnership. Vlahakis et al. I (12015> h 

SDP.81 was observed in October 2014 as part of 
ALMA’s Long Baseline Campaign, using between 22 and 
36 12 m-diameter antennas and ALMA’s band 4, 6 & 7 re¬ 
ceivers. The band 4 observations had the fewest total num¬ 
ber of antennas (a maximum of 27 compared to a maximum 
of 36 in the other bands), although the 21-23 element long 
baseline configuration was similar in all three bands. 

Four 1.875 GHz bandwidth spectral windows were used, 
over a total bandwidth of 7.5 GHz. In each observing band, 
one or two spectral windows covered a spectral line, with 
the remaining spectral windows used for continuum. The 
band 4, 6 and 7 data include the redshifted CO (5-4) (i> res t = 
576.267GHz), CO(8-7) (v rest = 921.799 GHz), and CO(10- 
9) (v r est = 1151.985 GHz) lines, respectively, as well as rest 

1 http://www.almascience.org 
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frame 250 /mi, 320 /mi and 500 /mi continuum. The band 6 
data also include the redshifted low-excitation water line 
H 2 0 (2o2 - 111) Orest = 987.927GHz; E up = 101 K) but 
we leave analysis of this feature for future work (see Section 

E}. 

The calibration and imaging of the data is described in 
the scripts provided on the ASP. These were carried out us¬ 
ing the C ommon Astronomy Software Application package 
fCASiFI: iMcMullin et al. |[20Q7h . A robust—1 weighting of the 
visibilities was used. Line-free channels were used to sub¬ 
tract the continuum emission from the CO data. The CO 
line data were imaged using rest frequencies correspond¬ 
ing to z — 3.042 and were uv- tapered to a resolution of 
~170mas (1000kA), since the high-resolution CO data has 
relatively low signal-to noise. In the case of H 2 0, the data 
were im-tapered to 200 kA (providing an angular resolution 
of ^0.9") in order to achieve a reasonable detection. The 
resulting RMS noise levels are 0.20 mJy and 0.15 mJy per 
21kms -1 channel for CO(5-4) and CO(8-7) respectively. 

We also carried out our own independent imaging of 
the calibrated visibility data supplied via the ASP. We at¬ 
tempted a variety of different tapers and cleaning parame¬ 
ters but found the ASP data to be already optimal for our 
purposes. We note also that in this paper we have used the 
band 4 image cube later staged on the ASP on March 2nd 
2015 with correctly subtracted continuum. 

For the purposes of our lens modelling, we binned the 
band 6 and band 7 continuum images from a pixel scale 
of 0.005" to a pixel scale of 0.01". This not only increases 
modelling efficiency, but also lessens image pixel covariance 
(see Section [3]) . In the modelling, we assumed the synthe¬ 
sised beam sizes prescribed in the ALMA data themselves; 
155x121 mas and 169x117 mas for CO(5-4) and CO(8-7), 
respectively. 

Panels (a) and (d) in Figure [I] show the ALMA band 6 
and band 7 continuum images. 

2.2 Near-IR data 

We have re-analysed the Hubble Space Telescope (HST) 
imaging of SDP.81 modelled by D14. We have applied our 
modelling to the deeper F160W image which has a total 
exposure time of 4418 s. The image was reduced using the 
IRAF MultiDrizzle package resulting in an image resam¬ 
pled to a pixel scale of 0.064". 

We have post-processed the data in two different ways. 
Firstly, we carried out an independent removal of the lens 
galaxy light prior to len s modelling using the GALFIT soft¬ 
ware (|Peng et al. 1120021 b In doing so, we have revealed ad¬ 
ditional structure in the F160W data to the south of the 
lens. Since this influences our new interpretation of the char¬ 
acteristics of the lensed source, we include this additional 
structure in our lens modelling with a larger image plane to 
encompass it (see section Ph2l for more details). 

Secondly, we applied a small astrometric shift of ~ 0.1" 
to align to the ALMA data. We determined this shift using 
the Sloan Digital Sky Survey data release 10 (|Ahn et al. I 

2 http://casa.nrao.edu 

3 The HST imaging was acquired with the Wide Field Camera 3 
(WFC3) in Cycle 18 under proposal 12194 (PI Negrello). 


l2014h to identify stars in the region covered by the HST 
image and then tied the r esulting stellar catalog ue to the 
two Micron All Sky Survey (|Skrutskie et al. 1120060 . We note 
that this astrometric alignment agrees very precisely with 
the independent alignment determined by shifting the best 
fit ALMA lens model centre to the centroid of the observed 
F160W flux. 

In the lens modelling of the F160W data, we used a 
point spread function created by the TinyTim software pack¬ 
age (jKrist I H993h . Panel (a) of Figure [2] shows the repro¬ 
cessed F160W data with the band 6 continuum overlaid as 
contours. 


3 LENS MODELLING 

We carried out our lens modelling in the image plane rather 
than the uv plane. There are two main advantages to this 
approach. The first is that the image can be masked to limit 
calculation of the goodness of fit to those parts of the image 
where there is detected emission from the lensed source. This 
gives a considerably more sensitive figure of merit for the 
fitting than working in the uv plane where modelling neces¬ 
sarily fits to visibilities that largely describe extended areas 
of background sky. The second is of particular relevance to 
the ALMA dataset under analysis in this paper; modelling 
in the image plane is vastly more efficient than working di¬ 
rectly with the extremely large visibility dataset which has 
to be trimmed in Fourier space anyw ay to en s ure t hat the 
modelling process is feasible fe.g. iRvbak et al.1l2015i ). 

The disadvantage of working in the image plane is that 
the pure interferometric visibilities are not directly mod¬ 
elled, but their modulated Fourier transform instead. A side 
effect of this is that image pixels are correlated by the beam 
which biases image plane modelling if the uncertainties do 
not take the covariance of image pixels into account. How¬ 
ever, in the case of the ALMA data modelled in this paper, 
the beam size is comparable to the image pixel scale and so 
this covariance is relatively low. The covariance is lowered 
even further by our use of the 2x2 pixel binned version of 
the band 6 and band 7 science verification continuum image 
data. Furthermore, the ALMA data have a very high cover¬ 
age of the uv plane which significantly reduces errors in the 
image plane. 

Any bias resulting from ignoring covariance in the im¬ 
age plane will affect the overall normalisation of the figure of 
merit. Whilst this effect will be small in the current ALMA 
data, this will still prevent a fair comparison between differ¬ 
ent lens model parameterisations in principle. However, for a 
fixed parameterisation such as that used in the present work, 
the relative difference in the figure of merit between differ¬ 
ent sets of parameter values remains unaffected. In this way, 
a reliable best fit lens model can still be found and hence 
image plane covariance is not a concern in this regard. 

We have verified that these assumptions are valid with 
the ALMA data analysed in this paper using the following 
procedure. Firstly, as we discuss below, we located the best 
fit lens model by application of the semi-linear method in the 
image plane to the 2x2 binned version of the band 7 contin¬ 
uum image. We then transformed the best fit model lensed 
image to the uv plane by using MIRIAD uvmodel to produce 
a simulated visibility dataset for the same uv coverage as in 
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Figure 1. Panels (a) to (c) show the reconstruction of the band 6 continuum image and panels (d) to (f) show the band 7 continuum 
reconstruction. The observed images are shown in panels (a) and (d), the model images of the reconstructed source are shown in panels 
(b) and (e) and the reconstructed sources are shown in panels (c) and (f). The white line in the source maps shows the position of the 
caustic. The inset zooms indicate the fidelity of the reconstruction. 


the ALMA dataset. Using the ALMA visibilities, their un¬ 
certainties and the model visibilities, we computed x 2 - We 
then varied different lens model parameters, stepping away 
from the set which provide the best fit in the image plane, 
generating new images each time and transforming to the 
uv plane to measure how % 2 varied. We found that although 
there is a slight offset in parameter space between the im¬ 
age plane minimum-x 2 and the uv plane minimum-x 2 , this 
is within the parameter uncertainties. 


3.1 Lens modelling procedure 

We used the late st implementation of the semi-linear in¬ 
version method (IWarren &; Dvel [2003) as described by 
Nighting ale fc Dve (1201.4) . The crux of the method is the 
manner in which the source plane discretisation adapts to 
the magnification produced by a given set of lens model 
parameters. By introducing a random element to the dis¬ 
cretisation and by ensuring that the discretisation maps ex¬ 
clusively back to only those areas in the image within the 
mask, the method eradicates biases in lens model parameter 
estimation. Crucially, the method removes a significant bias 


in the inferred value of the logarithmic slope of power-law 
mass density profiles which occurs when semi-linear inver¬ 
sion is used with a fixed sou rce plane size and/or a fixed 
source plane pixellisation fsee lNightingale &; Dve I [20151 . for 
more details). 

We have also used simultaneous reconstruction of multi¬ 
ple source planes from multiple images as described in D14. 
We attempted a variety of different combinations of data, 
but found that there were no clear improvements on lens 
model constraints beyond using a dual reconstruction of the 
ALMA band 6 and band 7 continuum image data. 

For our lens model, we adopted a single smooth power- 
law profile with a volume mass density of the form p oc 
r~ a . fn utilising this profile, it is assumed that the power- 
law slope, a , is scale invariant. This assumption appears 
to be reasonable on the scales probed by strong lensing as 
demonstrated by a lack of any trend in slope with the ratio of 
Einstein radius t o effective radius ([Koopmans et alTI 12006: 

iRuff et al. ]Sm|)- 

The corresponding projected mass density profile used 
to calculate lens deflection angles is therefore the elliptical 
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Figure 2. Source reconstruction of the HST F160W data. Panel (a) shows the observed image (note that the residual noise left from 
lens galaxy light subtraction is masked out in the lens modelling) with band 6 continuum contours overlaid and the inset panel shows 
the lensed image of the reconstructed source. In panel (b) we show the full scale of the reconstructed source that fits the tidal debris-like 
emission seen in the observed image. Finally, panel (c) shows a zoomed-in section of the source overlaid with band 6 continuum source 
contours. 


power-law profile introduced bv iKassiola &; Kovner I ( 199^ ) 
with a surface mass density, ft, given by 

ft = fto (r/lkpc) 1-a . (1) 

Here, fto is the normalisation surface mass density (the spe¬ 
cial case of a = 2 corresponds to the singular isothermal 
ellipsoid, SIE). The radius r is the elliptical radius defined 
by r 2 = x' 2 + y' 2 /e 2 where e is the lens elongation defined as 
the ratio of the semi-major to semi-minor axes. Three fur¬ 
ther parameters define the lens mass profile: the orientation 
of the semi-major axis measured in a counter-clockwise sense 
from north, 6 , and the coordinates of the centre of the lens 
in the image plane, (x c , y c )• Finally, following the findings of 
D14, we include an external shear component which is de¬ 
scribed by the shear strength 7 and orientation 0 1 measured 
counter-clockwise from north. 

As this paper is concerned primarily with the properties 
of the lensed source, we have not considered more compli¬ 
cated lens models. For example, one possibility is to model 
the lens using a cored density profile. Whilst there is contin¬ 
uum emission detected at the centre of the ring where a core 
would be expected to produce an image, our measurements 
of the spectral index indicate that this emission cannot be 
entirely from the background source. Similarly, the high res¬ 
olution and high signal-to-noise submm images may provide 
a perfect test-bed for lens mass profiles including substruc¬ 
ture. Nevertheless, we leave more detailed lens modelling for 
future study. 

4 RESULTS 
4.1 Lens model 

Table \l\ lists the most probable lens model parameters 
identified by marginalisin g over the Bayesian evidence (see 
iNightingale <V Dve 1120151 . for more details). The parameters 
are consistent with those obtained by D14 from modelling 
solely the HST data. They are also in agreement with the 


Lens parameter 

Value 

ftO 

(0.86 ± 0.04) x 10 10 M o kpc -2 

e 

1.25 ±0.04 

e 

13° ± 2° 

a 

2.01 ±0.05 

7 

0.04 ±0.01 

e 1 

-4° ± 3° 


Table 1. Most probable lens model parameters. The quantities 
are the lens mass normalisation, fto, the elongation of the lens 
mass profile, e, the orientation of the semi-major axis of the 
lens measured counter-clockwise from north, 0 , the density pro¬ 
file slope, a, the strength of the external shear component, 7 and 
the orientation of the shear 6 7 measured counter-clockwise from 
north. 

parameters obtained bv lRvbak et al. I (l2015h who performed 
uv modelling of the same ALMA data we have analysed in 
the present work. This provides further strength to our argu¬ 
ment concerning the viability of lens modelling in the image 
plane in this case. 

We also note that the ring is remarkably well fit if we 
force a slope of a = 2 and zero external shear, corresponding 
to a pure SIE model. In this case, the best fit elongation 
increases to e ~ 1.4 but the model has a lower evidence 
(with Ay 2 — 20) than the most probable model. 

4.2 Source morphology 

Figure |T| shows the lensed source reconstructed from the 
ALMA band 6 and band 7 continuum images. In the figure, 
we show the source reconstructed on a regular grid for pur¬ 
poses of illustration, although the adaptive source pixellisa- 
tion scheme was used in the lens modelling. We also show a 
zoomed-in section of the image of the source to demonstrate 
the quality of the fit. 

The source continuum emission shows that the distribu¬ 
tion of dust is very irregular with many small scale clumps. 
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Figure 3. Zeroth moment CO maps for C0(5-4) (panel a) and CO(8-7) (panel b) emission in the source, clipped at 2 cr. In both maps, 
the band 6 continuum emission is shown by the white contours and the F160W emission is shown by the yellow contours to illustrate 
the offset between the different wavelengths. 


The morphology is broadly similar between both bands, al¬ 
though the ratio of the band 7 flux to the band 6 flux varies 
slightly across t he so urce, mo st l ikely due to varying dust 
temperature. iRvbak et al I (|2015h also show that the star 
formation rate varies considerably across the source. 

Turning to the F160W data, Figure [2] shows our recon¬ 
struction performed using the best fit lens model obtained by 
simultaneously fitting to the band 6 and band 7 continuum 
data. As is immediately apparent from a direct comparison 
of the observed lensed images, there is a very distinct offset 
between the submm emission and the near-IR (rest frame 
blue) emission. The additional structure revealed by our re¬ 
processing of the F160W data is reconstructed in the source 
plane as shown in panels (b) and (c) of Figure [2j Panel 
(c) in the figure reveals the dominant near-IR component 
which gives rise to the bulk of the light seen in the ring. 
This matches that reconstructed by D14, although a tail of 
emission to the south is now apparent. In the larger scale 
reconstruction shown in panel (b), it appears that this tail is 
only part of a larger extent of near-IR emission. In addition, 
to the north of the main source component, there is another 
single component which, in the image plane, gives rise to 
the arclet seen to the north of the ring. We have measured 
the F110W-F160W colour of these additional components, 
drawing on the shallower F110W data described in N14, and 
find that it is consistent with the colour of the main arc seen 
in the HST data. 


We find total magnification factors of 16.0 zb 0.7 and 
15.8 =b 0.7 for the band 7 continuum and band 6 continuum 
reconstructed sources respectively. For the F160W source, 
we find a total magnification factor of 4.5 ±0.3 for the entire 
source as shown in panel (b) of Figure [2] For the dominant 
component seen in the F160W source, shown in panel (c) of 
Figure (2) we measure a magnification factor of 10.2 ± 0.5, 
consistent with the magnification measured by D14 for this 
part of the source. 


4.3 CO emission and source kinematics 

We used our best fit lens model to reconstruct the distribu¬ 
tion of flux in the source plane for each slice of the image 
data cubes released as part of the ALMA science verifica¬ 
tion data. Our reconstruction was able to detect significant 
CO(5-4) and CO(8-7) emission in the reconstructed band 
4 and 6 cubes respectively. We were unable to detect any 
significant CO (10-9) emission in the reconstructed band 7 
cube. 

Figure [3] shows the zeroth moment map of the CO line 
emission. These were made using CASA immoments, stack¬ 
ing over channels 31 to 61 inclusive (rest-frame velocities 
from -370 kms -1 to 260 kms -1 ) for CO(5-4) and channels 
30 to 60 inclusive (rest-frame velocities from -391 kms -1 to 
239 kms -1 ) for CO (8-7). The channel ranges were selected to 
fully encompass the spectral range of detected CO emission 
from the source. Our modelling yields total magnification 
factors of 13.7 ± 0.6 for the CO(5-4) flux and 13.1 dz 0.6 for 
the CO(8-7) flux. 

For comparison, in the zeroth moment maps we also 
show the the band 6 continuum reconstructed source (white 
contours) and the F160W source (yellow contours). It is im¬ 
mediately apparent from these plots that while the CO(5-4) 
emission follows the continuum emission, there is a signifi¬ 
cant offset between the more highly excited CO (8-7) emis¬ 
sion and the continuum. Furthermore, the CO(8-7) map 
shows more extended emission, linking the submm and near- 
IR emission regions. 

Using CASA immoments to generate a first moment map 
to obtain the velocity field in each reconstructed CO cube, 
we also found a relatively smooth variation in velocity across 
the source (see Figure 2]). The velocity field in each cube has 
the hallmarks of disk rotation and hence we fitted rotation 
curves assuming rotation of a disk (see Section I4.3.1D . Sim¬ 
ilar application of CASA immoments to generate the second 
moment map, shows a velocity dispersion which peaks in 
the dynamical centre of the source but also has strong peaks 
throughout. 
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4-3.1 Dynamical modelling 

To model the velocity field of the CO (8-7) and CO (5-4) and 
so estimate the disk inclination and asymptotic rotational 
velocity, we fitted a two dimensional model whose velocity 
field is described by a combination of stars and gas such that 
v 2 — v 2 D -\- v 2 h where the subscripts denote the disk and 
dark matter halo respectively (we ignore any contribution 
of Hi to the rotational velocity). For the disk, we assumed 
that the surface density follows a Fr eem an prof ile ( Freeman 
1970). For the halo, we assumed the lBerkert 1 ( 19951 ) density 
profile which incorporates a core of s ize ro and converges to 
the lNavarro. Frenk &; White I ( 1996 . NFW) profile at large 
radii. For suitable values of ro, the Berkert profile can mimic 
the NFW or an isothermal profile over the limited region of 
galaxy mapped by the rotation curves. 

This mass model has three free parameters: the disk 
mass, the core radius of the dark matter halo and the central 
core density. To fit the two-dimensional velocity fields, we 
constructed a two-dimensional kinematic model with these 
three parameters, but we also we allowed the [x / y\ centre 
of the disk, the position angle (PA) and the disk inclination 
to be additional free parameters. We constrained the [x / y\ 
dynamical centre to be within 0.5 kpc of the band 7 contin¬ 
uum centroid and then fit ted the two dimensional dynamics 
using an MCMC code fsee lSwinbank et al. 11201 2l . for further 
details) 

The best-fit kinematic maps and velocity residuals are 
shown in Figure [U The best fit disk inclination is (40 d= 
5)°. Assuming that the observed velocity field is indeed due 
to disk rotation, then the observed maximum line of sight 
rotational velocity of 210 kms -1 corresponds to an intrinsic 
asymptotic velocity of 320 kms -1 . 

Figure [5] shows the one dimensional rotation curves ex¬ 
tracted from a 0.4 kpc wide strip running along the major 
kinematic axis identified from the dynamical models for the 
CO(8-7) and CO(5-4) emission. For these rotation curves, 
we defined the velocity zero point using the dynamical cen¬ 
tre of the galaxy. The error bars for the velocities are derived 
from the formal lcr uncertainty in the velocity arising from 
the Gaussian profile fits to the CO emission. We note that 
the data have not been folded about the zero velocity, so 
that the degree of symmetry can be assessed. 

The rotation curves imply a mass within 1.5 kpc of 
(3.5 ±0.5) x 10 10 Mq with only a small contribution from the 
halo within this radius. Comparing the significantly lower 
maximum velocity dispersion of ~ 90 kms -1 with this ro¬ 
tational velocity suggests, under the assumption of a disk, 
that a correction for asymmetric drift need not be applied 
to the inferred dynamical mass. If we assume that all of the 
mass lies within the exponential disk component, this cor¬ 
responds to a surface mass density of 5000 ± 900M©pc -2 . 
T his is a typical valu e for high redshift SMGs (see Figure 6 
in llvison et al. 112013 ). 

Using this surface mass density and the observed ve¬ 
locity dispersion in the disk, the resulting Toomre stability 
value is Q ~ 0.3 which indicates a collapsing disk. Such 
low values of Q are observed in early merger systems b efore 
feedb ack can restore equilibrium in the disk (see [Hopkins] 

l2012h . 



Figure 5. Rotation curves of the lensed source from CO(5-4) 
and CO(8-7) emission, strongly resembling disk-like dynamics in 
both cases. The data points show the observed line of sight ve¬ 
locity extracted from a 0.4 kpc wide slit running along the major 
kinematic axis and the line shows the best dynamical model fit. 



Figure 6. Rest-frame de-magnified best fit SED to the band 
6 and 7 continuum source fluxes and to the fluxes from table 
2 in N14 (but de-magnified using our derived continuum source 
magnification factor). The dashed curves show the hot (99 K) and 
cold (48 K) SED components. 

4.4 Other intrinsic source properties 

4-4-1 Dust mass and total dust luminosity 

We have measured the dust mass of the lensed source by 
fitting a two-component spectral energy distribution (SED) 
model to a combination of our measured band 6 and band 7 
continuum fluxes and also the source fluxes presented in N14 
but de-magnified by our average continuum magnification 
factor of 15.9. The SED model allows the dust temperature, 
emissivity index, /3, and optical depth at 100 /xm, noo, to 
vary in the fitting. We used a dust mass opacity coefficient 
equivalent to ft 850 /xm = 0.077 m 2 kg -1 to be consistent with 
the fitting in N14. 

Figure [6] shows the best fit SED model which has an 
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Figure 4. Observed, modelled and residual velocity fields (panels (a), (b) and (c) respectively) and observed line of sight velocity 
dispersion (panel d) for the CO(8-7) line emission. The same quantities in the same order are shown for the CO(5-4) emission in panels 
(e) to (h). The white line in panels (a) and (e) shows the principle axis of the best-fit disk model and the dashed black line shows that 
of the stacked CO emission. The white ellipse in panels (b) and (f) shows the extent of the disk used in the dynamical modelling. 


effective dust temperature of 51 K, an emissivity index of 
/3 = 2.5 and an optical depth of noo = 6.5, yielding a dust 
mass of (1.8 ± 0.3) X 10 8 M© after de-magnifying by our 
average continuum magnification factor. The fit favours both 
a hot and cold dust component with temperatures of 99 K 
and 48 K respectively. 


If we fix the emissivity index to ft = 2.0, we obtain a 
lower optical depth of noo = 3 and a lower temperature of 
46 K but all optically thick fits, with or without fitting to 
the 160/xm flux (which corresponds to rest-frame ~ 40 /im 
and therefore not well fit by a modified black body SED) 
returned a dust mass in the range (1.8 — 2.1) x 10 8 M 0 . A 
second colder dust temperature component (10-40 K) was 
allowed in the fitting but was not favoured by the data. 


Assuming a typical gas to dust ratio of 150 (for exam¬ 
ple Dunne et al. 2000: [Praine et al~ 120071: ICoretese et al. I 


120121 : ISandstrom et al. 2013; S winbank et al. I [2014 ) there¬ 

fore gives a total molecular gas mass of (2.74=0.5) x 10 10 M 0 . 
In making these calculations, we have of course neglected 
any effects of differential magnification which could bias the 
inferred dust temperature and/or emissivity index. 


From integrating the best-fit SED between 8 and 
1000 /im, we estimate that the total far-infrared emission 
is (5.0 =L 0.8) x 10 13 L 0 , which agrees well with the value of 
5.4 x 10 13 L 0 given in N14. Our estimate of the average con¬ 
tinuum magnification factor of 15.9 implies that the intrinsic 
far-infrared luminosity is (3.1 d= 0.5) x 10 12 L 0 . 


44.2 CO(l-O) gas 

We have derived an estimate of the CO (1-0) luminosity of 
the gas in the lens ed source using the CO (1 -0) spectroscopic 
data acquired bv IValtchanov et al. I ( 2011 ). To do this, we 
firstly determined the spectrum of the total lensed flux (note 
that Valtchanov et al. plot the peak spectrum in their work). 
We then used our CO(5-4) source plane velocity and the 
magnification map from the lens model to estimate the mag¬ 
nification as a function of velocity in the source plane. By 
transforming the total CO (1-0) SED to velocity space, we 
then de-magnified the SED flux at each velocity with the 
corresponding magnification factor from our magnification- 
velocity relation. Figure [7| shows the de-magnified SED. 

Using the de-magnified SED, we determined a CO(1-0) 
flux of /co(i-o) = 85 =b 12mJykms -1 , which corresponds 
to a luminosity of L' co = (3.5 =L 0.6) x 10 10 Kkms -1 pc 2 . 
The errors take into account a number of uncertainties: 
1) There appears to be more CO (1-0) emission at higher 
and lower velocities than we see in our CO(5-4) data and 
so in our velocity-magnification relationship, we assumed 
a fixed magnification of 3 in these extremes; 2) The mea¬ 
sure of de-magnified flux is sensitive to the binning of the 
CO (1-0) spectrum and velocity-magnification mapping; 3) 
We have determined the de-magnified flux assuming the 
CO(1-0) emission exactly follows the CO(5-4) emission. The 
flux changes by ~ 15% if we assume the CO (1-0) emission 
is more evenly distributed in the source plane. 

To estimate the total mass of molecular gas in the 
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Figure 7. The de-magnified CO(1-0) spectrum in velocity space. 
The unshaded red histogram shows our estimated magnification 
factor. 

source, we used the C O luminosity to molecu lar gas mass 
conversion factor from iBothwell et al. I (|2013h of aco = 1 
(although there is a very large uncertainty on this factor; 
see for example IPapadopoulos et al. 1120121 ). This gives a to¬ 
tal molecular gas mass in the source of (3.5 ±0.6) x 10 10 M 0 
which is larger than the value of (2.7±0.5) x 10 10 M 0 we ob¬ 
tained by scaling from the dust mass, although both come 
with considerable uncertainty depending on the choice of 
assumptions (e.g., the dust-to-gas ratio, or the minimum 
amplification for gas beyond that seen in CO(5-4)). Never¬ 
theless, using the dynamical mass as an estimate of the total 
mass, and the range of gas masses derived from either the 
CO (1-0) or the dust luminosity, these values suggest a high 
molecular gas fraction in the central regions, in the range 
~ 75 — 100 per cent. 


5 DISCUSSION 

The offset between the reconstructed rest-frame optical 
source emission detected by HST/WFC3 and the recon¬ 
structed submm source from the ALMA data is strik¬ 
ing. D14 found good alignment between the rest-frame 
optical emission a nd th e submm source reconstructed by 
iBussmann et al.l (2013), although this latter study used 
imaging acquired with The Submillimeter Array (SMA) 
with a beam width approximately two orders of magni¬ 
tude larger than the ALMA image data analysed in the 
present work. The anti-correlation seen between the rest- 
frame optical and submm emission is not unique; the con¬ 
figuration in SDP.81 is similar to the offset between the opti¬ 
cal and submm emission seen in the z= 4.05 SMG GN20 (see 
iDaddi et al. f2 009: Ho dge et al. Il2015h . However, occurrence 
of such large offsets is not common. For exa mple, i n a sur¬ 
vey o f 126 SMGs carried out using ALMA (see lHodge et al. 1 
120131 ). only two sources, AL ESS88.il and AL ESS92.2 ex¬ 
hibit a similar configuration (|Chen et al. 1120151 ). 

There are a number of scenarios which could explain the 
observed configuration of the lensed source. The three most 
likely ones are: 1) The submm and optical components are 
simply two different sources which are closely aligned on the 


sky but well separated along the line of sight; 2) The submm 
and optical emission originates from the same source and 
we observe a total lack of optical emission from the submm 
region due to very strong dust absorption; 3) The optical and 
submm components are two separate systems undergoing a 
merger. 

The first scenario is challenging to verify, not least 
because of the large effective radius of the lens in the 
optical/near-infrared compared to the Einstein radius for 
the optical source. Photon noise from the lens light thus 
dominates the already faint optical source which will hinder 
attempts to measure its redshift, either spectroscopically or 
photometrically. This scenario also draws into question the 
structure of the optical source and the nature of the ap¬ 
parent tidal debris to the south and the component to the 
north. 

The second scenario could arise as a result of a strong 
dust lane in a disk galaxy. There are many examples of 
ULIRGs in the local Universe where strong absorption by 
dust lanes result in a delineation between optical and submm 
emission. In the case of SDP.81, it may be that a wide and 
thick dust lane is located at the eastern edge of the disk, 
which, because of its inclination, provides a much higher 
column density of dust toward regions in the source which 
emit optically. On the assumption that the optical absorp¬ 
tion efficiency factor of a dust grain is inversely proportional 
to wavelength, which is true in Mie theory if the grain size is 
much less than the wavelength, our estimate of the optical 
depth at 100 /xm (see Section 14. 4. ID implies that the opti¬ 
cal depth in the optical waveband is ~ 1000. Therefore, we 
would not expect to see much starlight from within the disk 
of gas and dust. 

The third scenario is suggested by the presence of the 
northern rest-frame optical component and the apparent 
tidal debris seen to the south. A tempting interpretation 
of this is that the northern component is a second galaxy 
which has passed through the larger, strongly lensed sys¬ 
tem, causing the tidal debris observed and enhancing the 
star formation rate and dust mass. Although the velocity 
field appears to be quite regular, the low Toomre Q param¬ 
eter we have measured (Q ~ 0.3) suggests a collapsing disk. 
Also, there are several strong peaks in the velocity disper¬ 
sion map which may still point towards some level of inter¬ 
action and there are filaments in the CO(8-7) emission (but 
not the less excited CO (5-4) emission) which extend up to 
the optical region. The strong emission in CO(8-7) is itself 
consistent with a merger scenario since this higher energy 
transition can not be produced over t he entire observed dis k 
by UV flux from star formation ([Papadopoulos et al. 11201 j ). 

In our attempts to search for more hints as to the na¬ 
ture of the source, we also carefully searched through all 
channels in the band 6 and 7 ALMA data cubes, looking for 
CO emission from the northern and southern optical com¬ 
ponents. Such emission might indicate an association with 
the source or provide additional kinematic information. We 
could not find any significant emission from the northern or 
southern optical components in either of the cubes. 

Whatever the connection between the rotating disk 
of gas and dust revealed by ALMA and the HST near- 
infrared sources, we can draw some conclusions about the 
properties of the star formation in the disk. We have es¬ 
timated the star-format ion rate in the disk from the in- 
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trinsic far- infrared l um inosity, which is the method sug¬ 
gested by iKennicutt <± Evans I ( 2012 ) for estimating the 
star-format ion rate in a highly obscured galaxy. From the 
relationship between star-format ion rate and far-infrared lu¬ 
minosity given in Kennicutt & Evans, we estimate that the 
star-format ion rate in this object is ~ 470 d= 80M©/yr. 

Although it is well-known that high-redshift galaxies are 
very c lump y an d irregu lar in broad-band optical/UV images 
(|Cowie. Hu fe Songalia Ifl 9951 ). it has always been an open 
question whether the clumpiness is the result of the star 
formation occurring in clumps or whether it is the result of 
patchy dust obscuration. In the case of this object, both the 
reconstructed CO and dust emission are clumpy on the scale 
of the point spread function in the reconstructed images, 
~ 200 pc. The CO lines are high-excitation lines, so that we 
cannot rule out the possibility that the clumpiness is the re¬ 
sult of a variation in the excitation rather than a variation in 
the distribution of the gas. There is also the recent sugges¬ 
tion that a clumpy CO distribution in high-redshift galax¬ 
ies might be th e result of th e destructio n of CO molecules 
by cosmic rays (iBisbas. Papadopoulos &; Viti I l2015h . How¬ 
ever, neither of these two caveats apply to the dust emission; 
dust grains are robust and found in all phases of the inter¬ 
stellar medium, and the emission from the dust depends 
only very weakly on the intensity of the interstellar radia¬ 
tion field. Therefore, the clumpiness of the dust is strong 
evidence that the distribution of gas in this object is truly 
extremely clumpy. The low value of the Toomre Q param¬ 
eter and the very irregular distribution of gas are exactly 
what one would expect if sections of the disk are collapsing 
to form stars. From the rotation curve and the velocity dis¬ 
persion, we estimate that the disk is unstable over the scale 
range ~ 50 pc to ~ 700 pc, the lower limit being the Jeans 
length and the upper limit being the scale on which the disk 
should be stabilised by shear. This agrees well with the sizes 
of the clumps observed. 

Finally, we have estimated the efficiency of the star- 
formation process in this galaxy. Using the dust mass as 
a tr acer of the to t al ma ss of gas, I Rowlands et al. I (|2014l ) 
and lSantini et al.~l (|2014l ) have found evidence that the star- 
formation efficiency (star-formation rate/gas mass) is higher 
in high-reds hift galaxies tha n in galaxies in the local Uni¬ 
verse. [Rowiands^r^r] ( 2014 ) give relationships between the 
star-formation rate and dust mass for local galaxies and for 
high-redshift SMGs. Using these relationships, we estimate 
that a typical galaxy in the local Universe with a dust mass 
equal to that of SDP.81 would be forming stars at a rate of 
7.2 ±1.3 M©/yr and that a typical SMG with the same dust 
mass would be forming stars at a rate of 58 ± 10 M©/year. 
Thus our estimated SFR of SDP.81 of 470M©/yr implies 
a star formation efficiency that is ~ 8 times greater than 
a typical SMG and ~ 65 times greater than in the nearby 
Universe. 


6 SUMMARY 

We have used the exceptional angular resolution of ALMA 
to reconstruct a detailed map of the submm emission and 
dynamics in the lensed source in SDP.81. Our modelling of 
the reprocessed HST data has revealed an offset of ~ 1.5 kpc 
between the submm and rest-frame optical centroids in the 


source. The submm continuum emission in the source is 
magnified by a total magnification factor of 15.9 ±0.7 which 
compares to the magnification of the rest-frame optical emis¬ 
sion of 10.2 ± 0.5 which mainly lies outside of the source 
plane caustic. Similarly, the CO(5-4) and CO(8-7) emission 
is magnified by the total magnification factors 13.7±0.6 and 
13.1 ± 0.6 respectively. 

Our reconstruction of the source kinematics from the 
CO emission reveals a relatively smooth velocity gradient 
across the source and suggests regular disk-like rotation. We 
have carried out dynamical modelling of the observed line 
of sight velocities and find that the data are best fit by a 
disk inclined at an angle of (40±5)° to the line of sight with 
an asymptotic rotational line of sight velocity of 210 kms -1 . 
Accounting for the disk inclination, this corresponds to an 
intrinsic asymptotic velocity of 320 kms -1 and an implied 
dynamical mass of (3.5 ±0.5) x 10 10 M© within 1.5 kpc. Our 
dynamical modelling returns a low Toomre Q-parameter of 
Q ~ 0.3. 

We have combined our measurements of the dust con¬ 
tinuum flux from the ALMA data with pho tometry of the 
lensed source given in iNegrello et al.l (|2014l ) to fit a mod¬ 
ified black body SED. This indicates a total dust mass of 
(1.8 ±0.3) x 10 8 M© after de-magnifying by our average con¬ 
tinuum magnification factor. Assuming a typical gas to dust 
ratio of 150 gives total molecular gas mass of (2.7 ± 0.5) x 
10 lo M©. We have also estimated the total molecular gas 
mass from t he de-magnffied GO(1-0) spectrum of the lensed 
source from I Valtchanov et al. I (|201ll ). This gives a total CO 
luminosity of L© 0 — (3.5 ± 0.6) x 10 10 Kkms -1 pc 2 which, 
assuming a gas mass conversion factor of unity, typical for 
ULIRGs in the local Universe, gives a total molecular gas 
mass of (3.5 ± 0.6) x 10 10 M©. 

One observable we have not discussed in this work 
is the low-excitation water line H 2 O ( 2 o 2 — I 11 ) (V rest — 
987.927 GHz ( E up — 101 K)) which can be seen in the band 
6 data. We have attempted to reconstruct the distribution 
of this line in the source plane using our most probable lens 
model, but this has proven too weak to locate easily. We 
have therefore left analysis of this feature for future study. 

To summarise, the nature of SDP.81 is somewhat per¬ 
plexing! Although the observational evidence we have assem¬ 
bled in this paper is suggestive of a galaxy merger, we can¬ 
not rule out other possibilities. Pending further analysis of 
additional observational data, the source evades our full un¬ 
derstanding. Nevertheless, this work has demonstrated the 
complexity we can begin to expect in high redshift SMGs 
when they are studied at the high angular resolution now 
made possible by ALMA’s incredible long baseline imaging 
capability. 
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Source property Value 


Total continuum magnification 
Total rest-frame optical magnification 
Total CO(5-4) magnification 
Total CO(8-7) magnification 
Disk inclination 
Asymptotic rotational velocity 
Toomre Q Parameter 
Total dust mass 

Total molecular gas mass from dust 
Dynamical mass within 1.5 kpc 
Total CO luminosity 
Total molecular gas mass from CO 
Star-formation rate 
CO(l-O) flux 

Total far-infrared luminosity 


15.9 ±0.7 
10.2 ±0.5 
13.7 ±0.6 
13.1 ±0.6 
(40 ± 5)° 

320 kms -1 
0.3 

(1.8 ±0.3) x 10 8 M© 

(2.7 ±0.5) x 10 10 M© 

(3.5 ±0.5) x 10 10 M© 

Leo = (3-5 ± 0 . 6 ) x 10 10 Kkms -1 pc : 
(3.5 ±0.6) x 10 10 M© 

470 ± 80 M© /yr 
34 ± 5 mJy kms -1 
(3.1 ±0.5) x 10 12 L© 


Table 2. Summary of the intrinsic (i.e. de-magnified where relevant) physical properties of the source. 
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